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Introduction 



Oil and water do not mix. Sometimes, however, we want them to mix and stay 
mixed for long periods of time, in the form of emulsions. The same applies to many 
other fluid pairs, A and B (say), including many that are molecularly miscible but only 
in some temperature range (usually high temperature). Quenching the A/B system 
initiates phase separation; these lectures describe that process and then various ways 
to halt it before it is complete, so that the resulting fluid domains have a finite length 
scale. 

Such finely divided mixtures, generally known as emulsions, have many applications 
in science and technology. Macro-emulsions, containing micron-sized or larger spherical 
droplets of oil in water or vice versa, are found in many products ranging from foods 
to agrochemicals. As we shall see these are generally not thermodynamically stable, 
but various tricks are available to prevent phase separation over long periods of time. 
Examples of such tricks include use of adsorbed surfactants to inhibit coalescence, or 
incorporation of trapped species to inhibit diffusive coarsening. 

For fluid pairs of similar phase volumes and viscosities, phase separation in 3D 
creates bicontinuous structure in which A-rich and B-rich phases each form connected 
interpenetrating domains. Stabilizing such a structure on a finite length-scale is more 
difficult than for droplets, although in principle this can be done by subjecting the 
system to continuous shear. A more practical alternative is to add a well chosen sur- 
factant, capable of reducing the interfacial tension effectively to zero. The resulting 
states, known as micro-emulsions, are thermodynamically stable and include not only 
droplet morphologies but bicontinuous ones. Bicontinuous micro-emulsions are fluids 
containing an interface whose fluctuating random geometry is sustained by entropy. 

Traditional surfactants adsorb and desorb easily and reversibly at the fluid-fluid 
interface, so that the interface in surfactant-stabilized micro- or macro-emulsions is 
generally in local equilibrium. This ceases to be the case for much larger amphiphilic 
objects, such as block copolymers, globular proteins, and janus particles, where in- 
terfacial detachment energies can be hundreds or thousands of times feT. Here the 
term 'janus beads' refers to spherical colloids with hemispheres of opposite surface 
chemistry; these obviously behave like large, irreversibly adsorbed surfactants. 

Perhaps less obviously, colloidal spheres of uniform surface chemistry can also be- 
come trapped, with similarly large attachment energies, at the fluid-fluid interface. 
Such structures have maximum local stability when the surface chemistry creates 
nearly equal interfacial tension between the solid and both fluids (neutral wetting). 
Stabilization by particles creates interesting alternative avenues to the formation of 
long-lived emulsified states. Such avenues were long neglected, possibly because to 



2 Introduction 



achieve them reproducibly requires detailed control of the entire preparation history 
of the sample. However this feature, which stems directly from the fact that detach- 
ment from the interface cannot be achieved by Brownian motion, is now seen as an 
opportunity for robustly locking in a desired microstructure. Such locked-in structures 
can include foams, multiple emulsions of nested droplets, and assemblies of highly non- 
spherical droplets. They also include bicontinuous states in which the jammed layer 
imparts rigidity to the entire 3D sample. 

In what follows, I will elaborate on the above narrative at the level of detail achiev- 
able in a set of four 90 minute lectures. This involves ruthless simplification, sometimes 
to the threshold of dishonesty (a necessity familiar to physics lecturers everywhere) 
which I hope I have not overstepped. Several calculations are left as exercises. (Also in 
some of these I have not carefully checked the numerical prefactors to the answers given 
- that forms part of the exercise!) Many important topics are deliberately avoided be- 
cause they are not central to the story: this applies particularly to systems where A 
and B are not simple fluids but themselves complex (polymers, suspensions, gels etc.). 

I have deliberately kept references through most of the text to the minimum level 
consistent with clarity. The exception is Section 5 on particle-stabilized emulsions 
which describes some relatively recent work. This area apart, there are good graduate 
texts and monographs that cover between them most of the topics addressed here at 
varying levels of detail; these include, but are not limited to, References [1; 2; 3; 4; 5; 
6; 7; 8]. 
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We start by recalling the description of an isothermal, incompressible, simple fluid 
with ncwtonian viscosity rj and density p. This obeys the Navier Stokes equation (NSE) 



Here the pressure field P must be chosen to enforce the incompressibility condition 



One generic approach to complex fluids is to consider a simple fluid obeying the NSE, 
coupled to a set of mesoscopic internal structural variables X(r), each a function 
of position. In principle the density and viscosity in (2.1) could depend directly on 
X(r). However a big simplification, which does not affect much the qualitative physics 
discussed below, is to assume these dependences are negligible. The main effect of 
X(r) is then to create an additional stress tensor that enters the NSE alongside the 
familiar viscous term: 



The stress term can alternatively be viewed as a force density f = Ver exerted by the 
structural degrees of freedom on the fluid continuum. 

This stress term can be large, and after suitable averaging its macroscopic effect 
in many complex fluids includes a greatly increased viscosity at small steady shear 
rates, and nonlinear rheological phenomena at larger ones. This applies for instance 
in both polymer solutions (where X describes some coarse-grained conformational 
variables for the polymers) and nematic liquid crystals (where X stands for the order 
parameter tensor Q a p describing local molecular alignment). For both these systems, a 
good numerical approach is to proceed via (2.3), supplemented by a recipe for <r[X(r)] 
and a time evolution equation for X itself. This approach means that a large body of 
numerical expertise in solving the forced NSE can be exploited. 

In these lectures we will be concerned with binary mixtures of simple Newtonian 
fluids, for which the relevant mesoscopic variable is a scalar describing the local com- 
position of the fluid mixture. We define it as 



p(v + v.Vv) = r/V 2 v - VP 



(2.1) 



V.v = 



(2.2) 



p(v + v.Vv) = ?7V 2 v - VP + V.er[X(r)] 



(2.3) 



0(r) = 



(n A - n B ) n: 



mcso 



(2.4) 



Here ua,b denotes the number of A,B molecules per unit volume locally and the over- 
bar denotes the macroscopic average over the whole system. The mesoscopic average 
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is taken over a large enough (but still small) local volume so that <p(r) is a smooth 
field. For simplicity we have assumed that A and B molecules have the same molecular 
volume; our assumption that p does not depend on <f> then requires that they also have 
equal mass. Given the incompressibility condition implicit in (2.2), we see that our 
composition variable obeys — 1 < < 1 with = 1 in a fluid of pure A. 
Later on we will show that for such a fluid mixture, 

SF 

V.o-(r) = -0(r)V/z(r) = "</>(r)V— -- (2.5) 

5<j)(r) 

where F[4>\ is a free energy functional (described further below) and its functional 
derivative /i is a chemical potential conjugate to the composition 4>. This is properly 
called the 'exchange chemical potential' as it controls the free energy increment on 
swapping B molecules for A (hence changing <j>) at fixed total density. 

2.1 The symmetric binary fluid 

We now restrict ourselves further to the case where the molecular A-A and B-B interac- 
tions are the same but there is an additional repulsion Eab between adjacent molecules 
of A and B. Combined with the previous assumptions about molecular mass and size, 
the system is now completely symmetric at a molecular level. At high temperatures, 
T > Tq — EAB/ks, the repulsive interactions are overcome by mixing entropy and 
the two fluids remain completely miscible. At lower temperatures however, the A-B 
repulsion will cause demixing into two phases, one rich in A, one rich in B. Entropy 
ensures that there is always a small amount of the other type of molecule present; 
close to the critical temperature Tc the two phases differ only slightly in <p. 

A schematic phase diagram for this system is shown in Fig. 2.1. The locus of co- 
existing compositions cf> = ±</>b(T) is called the binodal curve; within the binodal, the 
equilibrium state comprises two coexisting phases of composition ±0{,. The amount 
of each phase depends on the global composition (j) of the initial mixture. Specifically, 
the volumes occupied by the A-rich and B-rich phases, denoted $a,bV, where V is 
the overall volume of the system, obey §a + $b = 1 (clearly) and 

(<I>A - $ B )<t>b = 4> (2-6) 

Thus the phase volume $^4 evolves from zero to one as the overall composition <f> is 
swept across the miscibility gap from — to </>{,. Note that in what follows, we some- 
times refer to the coexisting fluids as simply A and B when the more proper terms 
would be A-rich phase and B-rich phase. Notationally however we distinguish by cap- 
ital letters phase volumes <&a,b (which lie between zero and unity) from compositions 
<j> (which obey — 1 < <f> < 1). 

The dotted line on the phase diagram is called the spinodal, <p — ±<fi s (T). Within 
this curve, a uniform initial state 4> is locally unstable to perturbations. On the other 
hand, between the spinodal and the binodal ((f> s < \<j>\ < <pb) the uniform state is 
metastable; to get started, phase separation requires nucleation of a large enough 
droplet. For a further discussion on the thermodynamics of binary fluid mixtures, see 

[i; 2]. 
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Fig. 2.1 Phase diagram of a symmetric binary fluid mixture. 



2.2 Mean field theory 

The simplest approach to the binary fluid is to postulate the following Landau- 
Ginzburg free energy functional 



where b and k are positive. Technically this is an expansion about a critical temperature 
Tc, at which the parameter a changes sign (positive above Tc, negative below). More 
importantly, apart from certain details at low temperature (where it does not predict 
correctly the exponential approach of 4>b — > 1), this functional can describe not only 
the physics of the phase diagram in Fig. 2.1, but also, in schematic form, the physics 
of interfacial tension. These issues can be tackled without further approximation by 
addressing (2.7) using field theory methods (the renormalization group being essential 
to understand the behaviour around Tc) but for our purposes, mean-field theory is 
sufficient. The mean-field theory is found simply by minimizing F . 

Two comments are in order. First, it would be possible to add a linear term J <j)dV 
to F. However, this equates to (ua — ^b)/(^a + n B ) — 4>V which is simply a constant 
governed by the global composition. Since it cannot vary, this term in F is physically 
irrelevant. Second, for a general fluid mixture one can expect a cubic term J (c^ 3 /3)dV . 
This term creates an asymmetric phase diagram, and is clearly important in fitting 
the model to real fluid pairs for which some asymmetry is always seen. However, for 
the physics discussed in these lectures, the cubic term adds lots of algebra and not 
much physics, so we omit it. 




(2.7) 
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The mean-field approach to (2.7) first considers states of uniform (f>(r) = <f>. For 
such states 

* = lP + \4> A = u{4>) (2.8) 

This has a single minimum at <f> — for a > 0, with positive curvature everywhere 
(Fig. 2. 2). The latter means that whatever <f> is chosen, one cannot lower the free 
energy by introducing a phase separation. On the other hand, for a < 0, F/V has 
negative curvature at the origin (and indeed everywhere between the spinodals cj) s = 
(—a/36) 1 / 2 ). Moreover it has two symmetric minima at <j> = ±<ftb with </>& = (—a/6) 1 / 2 . 
For \(f>\ < 4>b, F/V is minimized by demixing the uniform state at <f> into two coexisting 
states at <f> = ±<fib- (A small price must be paid to create an interface between these, 
but in the thermodynamic limit it is always worth paying.) The phase volumes of the 
bulk coexisting states is given by (2.6). 





Fig. 2.2 (a) Local part of free energy density U{(j>) for a > and a < 0. (b) The decay 
rate r(q) in the spinodal regime, (c) Resulting growth of a peak in the equal time density 
correlator S q (t). 



2.3 Interfacial tension 

These two bulk phases will organize themselves to minimize their mutual surface area; 
in most geometries, this requires the interface to be flat. To calculate its interfacial 
tension, we need to know the interfacial profile. We take a flat interface with its 
normal along the x direction so that 0(r) = (p(x). The boundary conditions are that 
<fr(x) approaches ±</> h at x = ±oo. To find the profile, we minimize F[tj>] — A / <j)dV with 
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these boundary conditions. (The A term ensures that the global composition remains 
fixed during the minimization.) The resulting condition 



S_ 



F — X <j)dV 



(2.9) 



shows that the chemical potential \i = SF/S(f> is equal to A and hence independent of 
position. Since for our symmetric choice of F\4>] we have /i = dU/d<fi — in a uniform 
bulk phase at density ±4>b, it follows that A = and 

fi(x) = acf) + be/) 3 - k\7 2 4> = (2.10) 

This expression for /i follows from the variational calculus [9], which evaluates the 
relevant functional derivatives as 
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^ J rdV = n^r)"" 1 (2.11) 
J (\7(f>) 2 dV = -2V 2 0(r) (2.12) 



It is then a mathematical exercise [1] to show that, with the boundary conditions 
already given, the solution of (2.10) is 

<j> {x) = ±<&,tanh(^^ (2.13) 



Co 

Here £ = (— «/2a) 1//2 is an interfacial width parameter, x marks the midpoint of the 
interface, and the overall sign choice depends on whether the A-rich or B-rich phase 
occupies the region at large positive x. 

Likewise one may show that the interfacial tension, defined as the excess free energy 
per unit area of a flat interface over the bulk phases, obeys [1] 

7o(z) = J ^ (y(M + U{<M) dx = {-^) (2.14) 

where 4> stands for 4>oi x ) obeying (2.13). The first two terms in the integrand give the 
full interfacial free energy evaluated for the equilibrium profile and the third subtracts 
off the bulk contributions. The interfacial profile is fixed by a trade-off between the 
penalty for sharp gradients (set by k) and the purely local free energy terms which, on 
their own, would be minimized by a spatial composition that jumps discontinuously 
from one bulk value to the other. All of the free energy parameters are ultimately 
controlled by molecular physics, but this dependence is system specific and we do not 
discuss it here. 

2.4 Stress tensor 



If the interfacial profile departs from the equilibrium one, a thermodynamic stress <r 
will act on the fluid. An important example is when the interface is not flat but curved; 
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under these conditions fj, cannot be zero everywhere. For use in the NSE we require 
not the stress tensor directly but the thermodynamic force density f = V.cr which, as 
stated in (2.5) can also be expressed as —tpVfx. 

One derivation of this result is to consider moving a small blob of fluid of volume 
AV and composition (f>, from one part of the system to another. Denoting the chemical 
potential at its original position by /i and that at its final position by /i + Sfx, the free 
energy change SF is (arguably) <j>5/j,AV. If this happens within a region where fj, varies 
slowly so that 5fi = V/J,.5r, we have 

SF = AV(cf)Vn).6r (2.15) 

This is the work done by an external agent who must therefore apply a force density 
Thus the internal force density exerted by the surroundings on the blob is 

f = -0V^ (2.16) 

which is the required result. 

The above argument is seductive but incomplete: we are working at fixed fluid den- 
sity, so some unspecified other blob of fluid must move in the opposite direction from 
the one we have transported. A more correct, but much longer, derivation is to con- 
sider an incompressible deformation of the binary fluid mixture, in which line elements 
transform as dri — > dr.i + eijdrj, with the strain tensor. This deformation moves 
fluid elements about but cannot change their local composition: if a material point 
moves from r to r' then 4>(r') after deformation equals <fi(r) before. This information 
is sufficient to allow the free energy functionals F[<j>] before and after the deformation 
to be compared; to leading order in one has 

AF=^a ij e ij (2.17) 
The resulting stress tensor is found (eventually) to be [10] 

Sij - K (di4>)(d j( p) (2.18) 

from which (2.16) can be confirmed (this is left as an exercise). 

2.5 Equation of motion for composition; Model H 

Having specified how to compute the force term in the NSE for a binary fluid from 
the composition field <ft(r), we now need an equation of motion for this quantity itself. 
This takes the form 

(f> + v.V</> = -V.J (2.19) 

where the term in v represents advection by the fluid velocity, so the left-hand side is 
the co-moving time derivative of (f>. This derivative must be the divergence of a flux, 
because A and B particles are not created or destroyed and thus cf) is a conserved field. 
The form chosen for the flux is 

J = -M(<j>)Vn (2.20) 

where M(4>) is a mobility that depends locally on composition. More generally one 
could have a nonlocal M [</>]; but for our purposes it will generally be enough to consider 



^0 4 + K 0V 2 0+^(V0) 
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M{<p) to be constant. This mobility factor describes, under conditions of fixed total 
particle density, how fast A and B molecules can move down their respective chemical 
potential gradients to relax the composition field. 

Combining (2.19,2.20) with our earlier results for the chemical potential and the 
forced NSE, we arrive at a closed set of dynamical equations with which to describe the 
dynamics of phase separation in an incompressible, isothermal, binary fluid mixture: 



H{r) = a<j)(r) + 6</> 3 (r) - nV 2 (f>(r) (2.21) 

p(v + v.Vv) = ryV 2 v - VP - cjN y, (2.22) 

V.v = (2.23) 

4> + v.V<t> = V.(MVp) (2.24) 



Collectively the above equations are known as 'Model H' [1]. Note that the pressure 
field P has the job of ensuring incompressibility so that (2.23) is obeyed. In this 
constraint-enforcing role, P takes on whatever value it is told to by the other terms 
in the equations and cannot therefore introduce any new scaling behaviour that was 
not already calculable from those terms. 

As we have derived them, the equations of Model H are of mean-field form: they 
are deterministic, and take no account of noise. Noise terms are however important in 
at least two situations. One is near the critical point (not addressed in these lectures) 
where thermal fluctuations play a dominant role in the statistics of <\>: the mean- 
field theory implicit in the noise-free treatment breaks down. Another case where 
noise matters is when one has suspended fluid droplets (and/or colloidal particles). 
These objects will move by Brownian motion, which arises from thermal momentum 
fluctuations. These are neglected by our noise-free NSE, so that droplets of one fluid 
in another cannot diffuse. 

Fortunately the noise terms can easily be determined using the fluctuation dissi- 
pation theorem [1]. For the order parameter fluctuations we need to add to (2.24) a 
term —V.J™ where the random current J™ (superscript n for noise) has the following 
statistics: 

(J?(r,t)J?(r',t')) = 2k B TM5 ij S{r-r')S(t-t') (2.25) 

Similarly, to include Brownian motion we need to add to (2.22) a term — V.er™ where 
er™ is a fluctuating thermal stress whose statistics obey [11] 

(^■MK^iV)) - 2k B Tr,6 ik 5 jl 8(r - r')5(t - t') (2.26) 

With these terms added, Model H is transformed from a mean-field approximation to 
a complete description of the binary fluid whose free energy functional is F[4>]- Accord- 
ingly, in the absence of flow and with noise terms included, Model H ultimately achieves 
in steady state the probability density V[<j)\ oc exp (—/3F [(/>]), with f3 = (fe^T) -1 , as 
required by the Boltzmann distribution. 
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Having assembled the conceptual and mathematical tools we need, let us now examine 
some of the dynamics of phase separation. 

3.1 Spinodal decomposition 

We start by considering the spinodal instability. Ignoring advection initially, we write 
(2.24) as 

^V.(MV^i) (3.1) 
= V. (MV [U'{4>) - kV 2 0] ) (3.2) 
= V. (M [U"((f>) V<j> - kV 2 0] ) (3.3) 

Here £/(</>) is the local free energy density of a uniform state as defined in (2.8), and 
primes denote differentiation of this function with respect to <p. 

Next we linearize this equation about a uniform initial composition <f>, and fourier 
transform, to give 

<j> q = -M(4>)q 2 [U"(4>) + nq 2 ] q> q = -r{q)4> q (3.4) 

where the second equality defines a wave- vector dependent decay rate r(q). For U"(<f>) > 
this is positive for all q: all fourier modes decay and the initial state is stable. In 
contrast for U"((j)) < 0, the system is unstable, with r(q) negative at small and inter- 
mediate wavevectors. (Stability is restored at high enough q by the n term.) Differ- 
entiating the growth rate — r(q) with respect to q we can identify the fastest growing 
instability to be at q* = -U"{4>)/2k (Fig.2.2). 

Even neglecting the noise terms in the dynamics (2.25,2.26), the initial condition 
can be assumed to have some fluctuations. Those whose wavenumber lies near q* grow 
exponentially faster than the rest, so that the time dependent composition correlator 
S q (t) = (4> q (t)(j)- q (t)) soon develops a peak of height scaling as exp[|r(g*)|i] around 
q* (Fig.2.2). Hence during this 'early stage' of spinodal decomposition a local domain 
morphology is created by compositional diffusion (inter-diffusion of A and B) with 
a well defined initial length scale set by ir/q* . The amplitude of these compositional 
fluctuations grows until local values approach ±</>&, the composition at which A-rich 
and B-rich phases can coexist. There soon develops a domain pattern, still initially 
with the same length scale, consisting locally of these phases, separated by sharp 
interfaces whose local profiles resemble (2.13). 
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3.2 Laplace pressure of curved interfaces 

As mentioned previously, however, unless these interfaces are perfectly flat, they will 
exert forces on the fluid via the —cfiVp: term in (2.22), in response to which the fluid will 
be set in motion. The physics of this term, for interfaces that are locally equilibrated 
but not flat, is that of Laplace pressure. 

To remind ourselves of this physics, consider a spherical droplet of one fluid in 
another with radius R and interfacial tension 7. Let the pressure inside the droplet be 
greater than that outside by an amount AP. The total force on the upper half of the 
droplet exerted by the bottom half is then 



which must vanish if the droplet is not moving. The first term comes from the vertical 
component of the extra pressure acting across the equatorial disc, and the second is 
the tension acting across its perimeter. Hence AP = 2j/R; equilibrium requires the 
internal pressure to be higher as a result of curvature and interfacial tension. More 
generally one has a Laplace pressure 



where R\ and i? 2 (see Section 4.4) are the principle radii of curvature of the interface. 

3.3 What happens next? 

The next stage of the phase separation kinetics depends crucially on the topology of the 
newly formed fluid domains. This is controlled mainly by the phase volumes $a,b of the 
A-rich and B-rich phases. Roughly speaking, if 0.3 < $>a < 0.7 the domain structure 
will be bicontinuous: one can trace a path through the A-rich phase from one side of 
the sample to the other, and likewise for the B-rich phase. Outside this window, the 
structure instead has droplets of A in B ($,4 < 0.3) or B in A ($^4 > 0.7). The values 
0.3 and 0.7 are rule-of-thumb figures only, with details depending on many other factors 
(including any asymmetry in viscosity [2; 12]) that we do not consider here. Note also 
that the window of bicontinuity shrinks to a single point in two dimensions, where the 
slightest asymmetry in phase volume and/or material properties will generally result 
in a droplet geometry in which only one phase is continuous. 

3.4 Coalescence of droplet states 

If the post-spinodal structure is that of droplets, each relaxes rapidly to minimize its 
area at fixed volume resulting in a spherical shape. After this, Laplace pressures are 
locally in balance and although stresses are still present there is no net fluid motion 
in the absence of noise. Thermal noise however allows droplets to explore space by 
Brownian motion, and the resulting collisions cause the mean droplet radius R to 
increase by coalescence. To estimate the rate of this process, we observe that for A 

— 1/3 

droplets in B, the mean inter-droplet distance L is of order R$ A at small $,4 or, 
more generally, L ~ Rf{<f>A)- Each droplet will collide with another in a time At 



nR 2 AP - 2tt 7j R = 



(3.5) 




(3.6) 
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of order L 2 /D where D ~ ksT/riR is the diffusivity. Upon collision, two droplets of 
radius R make a new one of radius 2 1 ^ 3 R causing an increment Ami? = (ln2)/3. This 
gives 

Ami? k B T 
At ?7i? 3 v ; 

where the left hand side can now be approximated as dhiR/dt — R/R. By integration 
we then obtain the scaling law 

«<o~(^r 0.8) 



This argument assumes that coalescence is diffusion-limited, and shows that in this case 
Brownian motion will cause indefinite growth of the mean droplet size, culminating in 
total phase separation. 

The assumption of diffusing spherical droplets is reliable at low phase volumes 
$,4 of the dispersed phase, but when this is not small more complicated routes to 
coalescence, some involving droplet-scale or macroscopic fluid flow, are possible. One 
of these is so-called 'coalescence-induced coalescence' where the shape relaxation post- 
collision of a pair of droplets creates enough flow to cause another coalescence nearby 
[13]. This gives a new scaling (R <~ "ft/r/) which coincides with one of the regimes 
described later for the coarsening of bicontinuous structures ((3.19) below) and indeed 
stems from the same balance of forces as will be discussed for that case. 

In many droplet emulsions it is possible to inhibit the coalescence step, so that this 
route to phase separation is effectively blocked. For instance, adding charged surfac- 
tants can stabilize oil droplets in water against coalescence by creating a coulombic 
barrier opposing the close approach of droplet surfaces. Steric interactions between 
surfactant tails can likewise stabilize water-in-oil emulsions. If the rupture of a thin 
film of the continuous phase between droplets has a high enough nucleation barrier, 
coalescence rates can often be reduced to a manageable or indeed negligible level [3] . 



3.5 Ostwald ripening 

Sadly however, switching off coalescence is not enough to prevent macroscopic phase 
separation of droplet emulsions. This is because of a process called Ostwald ripening, in 
which material is transported from small droplets to large ones by molecular diffusion 
across the intervening continuous phase. 

This process is most easily considered when (fif, is close to unity so that A droplets 
consist of nearly pure A. There is nonetheless a small equilibrium concentration c eq 
of A molecules in the nearly pure B phase: this is what allows diffusion of A between 
droplets. 

The driving force for Ostwald ripening is the Laplace pressure difference between 
small and large droplets. In a droplet of radius R the Laplace pressure H = 2-f/R causes 
the chemical potential ha of A molecules within it to be raised by A/j,a(R) = TIva 
where va is a molecular volume. Because the interface is in local equilibrium, the 
chemical potential just outside a droplet of radius R is also raised by this amount. 
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Treating the dilute solution of A in B exterior to the droplet as an ideal mixture, we 
then have a concentration field ca{t) at radius r from the droplet centre that obeys 

CA(R + ) = c eq {l + ^) (3.9) 

This concentration is high outside small droplets and lower outsider large ones; the 
resulting gradient causes a diffusive flux from smaller to larger droplets. 

To address this at mean field level we consider just one droplet, of radius R, obeying 
(3.9) and take a boundary condition at infinity 

c A {r -> oo) =c(t) =c eq (l + e(t)) (3.10) 

Here the e(t) term denotes the fact that the system globally has not yet reached 
equilibrium: hence there is a mean supersaturation of A in B resulting from the fact 
that distant droplets are themselves of finite radius. (Note that e(t) will tend to zero 
if and when the mean droplet size tends to infinity.) 

It is a simple exercise then to solve the quasi-steady diffusion equation DaV 2 ca = 
with the boundary conditions (3.9,3.10) to find 

c A (r) = c(t) + R( C{R+) ~ 5 ) (3-11) 



and from this to derive the flux of A molecules onto (or off of) the droplet surface. 
The result is an equation for the droplet size 

v A D A c eq ( 2-fv A 



^ q (312) 

The function R(R) is shown in Fig. 3.1. This exhibits an unstable fixed point at 

= (3-13) 

Droplets bigger than this grow, those smaller, shrink. 

We can now find the scaling of the typical droplet size by assuming this to be 
comparable (but not exactly equal) to R e : 



R( t ) ^ c^y )" 3 ~ t i " (3.i5) 



giving the scaling law [3] 



k B T 

From this it follows that the global supersaturation varies as e(t) <~ t -1 / 3 . A more com- 
plete theory, due to Lifshitz and Slyozov, not only confirms these scalings but gives 
detailed information on the droplet size distribution [14]. Note that (3.15) has similar 
time dependence to (3.8) for coalescence; this stems from the fact that both mecha- 
nisms arc ultimately diffusive. However, the nature of the diffusing species (droplet in 
one case, molecule in the other) is quite different, resulting in prefactors that involve 
unrelated material properties for the two mechanisms. 
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Fig. 3.1 Growth rate of droplet as a function of size during the Ostwald process (a) without 
and (b) with trapped species. 

3.6 Preventing Ostwald ripening 

We see from (3.15) that the Ostwald process can be slowed by reducing the interfacial 
tension 7. This is discussed later in relation to surfactants, but unless the tension 
is reduced effectively to zero this will only slow things down by a moderate factor. 
One might also reduce the solubility c eq of A in B, but generally the aim is to make 
emulsions of particular fluids, so replacing A with a less soluble species is often an 
impractical suggestion. On the other hand, a closely related approach is to include 
within the A phase a modest concentration of a species that is effectively insoluble 
in B. This might be a polymer or, if A is water and B oil, a simple salt. The idea is 
that the trapped species in the A droplets creates an osmotic pressure which rises as 
R falls, hence opposing the Laplace pressure. Treating the trapped species as an ideal 
solution in A, (3.12) is replaced by [15] 

• v A D A c eq ( 2jv A (v A \ 

R= ^R—{ € -^fR + WMw) (3 ' 16) 

where the last term is Tlov A /kBT with IIq the osmotic pressure of £ trapped particles 
within a droplet of radius R. Figure 3.1 shows the new structure of the R(R) equation. 
There is now a stable fixed point at a size determined by the trapped species 




(3.17) 



where i?o is the initial droplet size and ct the initial concentration of trapped species in 
the A phase. Droplets that have shrunk to a size R^(Ro) can coexist with a bulk phase 
of A without shrinking further: the Laplace pressure is balance by n<> Indeed if the 
initial size obeys i? < R((Ro) (treating the droplets as monodisperse for simplicity) 
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the Ostwald process is switched off entireiy. By this route one can thus make robust 
'mini-emulsions' [16] or 'nano-emulsions' [17], which will not undergo coarsening by 
the Ostwald process. However these are still metastable: so long as the tension 7 is 
positive, the free energy can always be reduced by coalescing droplets to reduce the 
interfacial area. 

3.7 Coarsening of bicontinuous states 

As mentioned previously, for 0.3 < $a < 0.7 (roughly speaking) the domains of A-rich 
and B-rich coexisting fluids remain bicontinuous. This allows coarsening by a process 
faster than either coalescence or Ostwald ripening, in which the Laplace pressure 
gradients pump fluid from one place to another. The driving force is interfacial tension 
and at any time there is a characteristic domain length scale L(t), much larger than the 
interfacial width, which we shall assume to be the only relevant length in the problem, 
so that V <~ 1/L. Moreover the characteristic magnitude of the Laplace pressure is 
n <~ j/L(t). This fixes the scale of the forcing term in the Navier Stokes equation 
(2.22) as — (j>Vfi ~ VII <~ j/L 2 . The fluid velocity is of order L so the viscous term 
scales as r/L/L 2 . The inertial terms are pir ~ pL and pv.Vv ~ p(L) 2 /L. The VP 
term, which ensures incompressibility, is slave to the other terms (see Section 2.5). 

An important property of (2.22), once sharp interfaces are present so that </>V/j ~ 
7/L 2 , is that it contains only three parameters, p, 7,77. From these three quantities 
one can make only one length, L = ?? 2 /p7, and one time to — i] 3 /pj 2 . This means 
that the domain scale L(t) must obey [18; 19; 20] 

m = f (t_\ 

L T \toJ 

where, for given phase volumes, f(x) is a function common to all symmetric binary 
fluid pairs. Note that noise is excluded, so that one expects this scaling to fail in droplet 
regimes where diffusion is important. (Noise may also be important in 2D [21].) It can 
also fail in bicontinuous states at relatively early times and very high viscosities when 
an Ostwald-like process can dominate over fluid flow, giving i 1 / 3 scaling as in (3.15). 

Excluding that regime, in bicontinuous states we therefore expect different be- 
haviour according to whether L/L is large or small. For L/L small it is simple to 
confirm that the inertial terms in (2.22) are negligible. (Note also that, in any regime 
where f(x) is a power law, the two inertial terms have the same scaling.) The primary 
balance in the NSE is then r]L/L 2 <~ 7/L 2 resulting in the scaling law L(t) <~ jt/t] so 
that 

f(x) (xx ; x < x* (3.19) 

This is called the viscous hydrodynamic or VH regime [18]. In contrast, at large x the 
primary balance in the NSE is between the interfacial and inertial terms. It is a simple 
exercise then to show that L(t) <~ ^/p) 1 / 3 * 2 / 3 so that 

f(x) cx x 2/3 ; x > x* (3.20) 

This is called the inertial hydrodynamic (IH) regime [19]. In (3.19,3.20) we have in- 
troduced a crossover value x — x* between the VH and IH regimes. In practice this 



(3.18) 
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crossover is very broad, and the crossover value rather high: x* ~ 10 4 . The high 
crossover point is less surprising if one calculates a domain-scale Reynolds number 

H.-f -/(.)£ (3.2D 

The crossover value of Re then turns out to be of order 10 [20], and the largeness of 
x* is found to stem from a small constant of proportionality in (3.19). It means that 
in practice a clean observation of the IH regime has only been achieved in computer 
simulation: in terrestrial laboratory experiments the domains are by then so large 
that the slightest density difference between A and B causes gravitational terms to 
dominate. 



3.8 Shearing binary fluids 

The coarsening of bicontinuous demixed states described above leads inexorably to 
complete phase separation; in practice this is something that, as explained in the 
introduction, we often wish to avoid. In a processing context, it is sometimes enough 
to temporarily maintain a well-mixed, emulsified state merely by stirring the system. 
Though industrial stirring is complicated, for our purposes it is enough to consider 
the effects of a simple shear flow. 

Consider such a flow with macroscopic velocity along x, and its gradient along 
y; z is then the neutral (or vorticity) direction. In simulations one can use boundary 
conditions with one static wall at y = 0, another sliding one at y = A and periodic BCs 
in x, z - or in practice there are ways to introduce periodic BCs also in y. Nonetheless, 
the system size in that direction, A is important in what follows. The top plate moves 
with speed A/t s where l/t s (usually denoted 7) is the shear rate. 

The question we ask is whether nonequilibrium steady states now exist for which 
the fluid domains have finite length scales L x , VtZ in all three directions. The simplest 
hypothesis is that these lengths, if they exist, all have similar scaling: L XtVtZ ~ L. 
Moreover, given the preceding discussion of terms in the NSE we expect in steady 
state that L/L is now a function, not of t/t but of t s /t . That is, the previous 
dependence on time is no longer present (because we assume a steady state exists) 
but is replaced by a dependence on the inverse shear rate. The functional form of this 
dependence could in principle be anything at all, but the simplest scaling ansatz is 
that the system coarsens as usual until t ~ t s , whereupon the shearing takes over and 
L stops increasing. If so, 

(3 - 22) 

where f(x) is the same function as introduced previously, for which (3.19,3.20) hold. If 
so, for t s /to <C x* we have L/L ~ t s /t a and for t s /t » x* we have L/L ~ (t s /t a ) 2 ^ 3 . 
The Reynolds number obeys Re^ f(x)df/dx as in (3.21). In contrast to what happens 
for any problem involving shear flow around objects of fixed geometry, Re is now small 
when the shear rate is large and vice versa. This is because at small shear rates, very 
large domains are formed. 
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The above picture is the simplest possible [22] . We note only two of several possible 
complications. First, in principle L x y ^ z could all have different scalings. The resulting 
anisotropies could spoil any clear separation of the VH and IH regimes; with three- 
way force balance in the NSE, there is no reason to expect clean power laws for any of 
these quantities. Secondly, at high shear rates the system-size Reynolds number ReA ~ 
pA 2 /r]t s becomes large. The presence of a complex microstructure could promote any 
transition to conventional fluid turbulence expected in this regime. 

In practice, experimental tests of the above predictions are patchy (see [2] and 
references cited in [23]). Gravity complicates matters, as does viscosity asymmetry 
(unavoidable in practice) between phases. However, relatively clean tests are possible 
via computer simulation [23; 24]. In 2D one finds apparent scaling laws L x /L ~ 
(t s /to) 2 / 3 and Ly/Lo <~ (t s /to) 3 / 4 over a fairly wide range of length and timescales, 
most of which are however in the crossover region around x* [24] . The fitted exponents 
change slightly if instead of the flow and gradient direction one uses the principal axes 
of the distorted density patterns, but still do not coincide. In 3D, where the simulations 
require very large computations, the IH (2/3 power) scaling has been observed within 
numerical error for all three length scales within a range of accessible domain-scale 
Reynolds numbers Re^ between 200 and 2000 [23]. These measurements are however 
limited by the onset of a macroscopic instability to turbulent mixing at ReA — 20, 000. 

A snapshot of the highly distorted domain structure seen in computer simulations 
of sheared binary fluids (within the laminar flow regime) is shown in Fig. 3.2. 




Fig. 3.2 Image of binary fluid domains in a nonequilibrium steady state created by shearing. 
This is a 2D simulation, but a slice through a 3D run looks very similar. See [23; 24] . 



4 

Stabilizing Emulsions 
Thermodynamically 



Surfactants, comprising small amphiphilic models with a polar head group and an 
apolar tail, are well known to reduce the interfacial tension between coexisting phases 
of oil and water (as well as other pairs of apolar and polar fluids). Moreover they 
generally have fast exchange kinetics between the interface and at least one bulk phase 
in which they are soluble; this means that the interface remains locally in equilibrium, 
even if at larger scales the emulsion is only metastable. Use of surfactants can thus be 
viewed as a thermodynamic route to the stabilization of interfacial structures. 

4.1 Interfacial tension in the presence of surfactant 

The basic effect can be illustrated by first considering an ideal solution of surfacant 
molecules each carrying a unit polarization vector p 4 . (This simply denotes orientation; 
a coulombic dipole is not required.) A local coarse graining creates a smooth field 
p(r) = (pi) meso ; because the solution is ideal, it is easy to show that the variance 
X = (|p| 2 ) 01 the fluctuating p field is proportional to the global concentration c s of 
the surfactant molecules. Therefore for noninteracting (ideal solution) surfactants one 
can write down a free energy to describe these fluctuations as 

where (in terms of a concentration-independent parameter %) 

X(cs) = xc s (4.2) 

is an osmotic compressibility. We next add a coupling term to represent the reduction 
in free energy caused when a surfactant molecule resides at the A-B interface with its 
polarity suitably aligned along the composition gradient V</>: 

F c = J vp.V<j)dV (4.3) 

The remaining terms of the free energy are simply those of the Model H binary fluid 
mixture, from (2.7), so that we now have 

F[4>,p]= J ^(^ 2 + ^ 4 + ^(V0) 2 + ^|p| 2 + ^p.V^ (4.4) 



(4.1) 
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It is a simple exercise to minimize this over p(r) at fixed 4>(r), and a slightly more com- 
plicated one to explicitly integrate over the fluctuating p field by Gaussian integration 
to obtain e~P F ^ = J e~^ F ^' p ^T>p. The result of either calculation is to recover the 
original Model H free energy (2.7), but with a renormalized square gradient coefficient 



(4.5) 



From this it follows that the interfacial tension varies with the concentration of an 
ideal surfactant as 



7(c s ) 



-8a 3 (k — v xc s ) 
%2 



1/2 



(4.6) 



This vanishes, with infinite slope, at a concentration c s = c = k/p 2 x (Fig. 4.1). 

* Y 




C* c 



Fig. 4.1 Interfacial tension as a function of surfactant concentration: (a) for an ideal surfac- 
tant solution, and the case where c* > c; (b) for a typical surfactant where c* < c, so that 
micellization pre-empts the vanishing of the interfacial tension 7. 

The above calculation for an ideal solution of surfactant gives a simplified picture 
of how, on increasing surfactant concentration, the interfacial tension between two 
immiscible fluids A and B can vanish. (For a far more advanced discussion, see [25].) 
Before considering what happens when this point is reached, we must discuss why, in 
practice, it very often is not reached. This is because surfactant solutions are nonidcal, 
due to the phenomenon of micellization in which individual molecules become aggre- 
gated into micelles containing several tens of molecules. The effect of this process, as 
we shall see below, is to put a cap on the osmotic compressibility x( c s)- This happens 
at a certain concentration c* ; only if this lies beyond c as defined above does 7 come 
close to zero. 



20 Stabilizing Emulsions Thermodynamically 



4.2 Micellization 

As just explained, surfactant solutions are not ideal. First recall that in an ideal 
solution, the concentration and chemical potential are linked by 

c s = v^ 1 exp[/3/z s ] (4.7) 

where vt is a molecular volume. (This is the same expression as for an ideal gas where 
vt is instead the thermal de Broglie volume.) Next, note that the role of x is to tell us 
the cost of pulling surfactant molecules out of solution to put them on the interface: 
this is indeed the meaning of (4.1). It should be clear therefore that for a nonideal 
solution what matters is not c s but fj, s : we therefore expect (4.1) to be replaced by 

X(c) =Xv T le MPHs] (4.8) 

The remaining task is to understand how fi s behaves in the presence of micellization. 

The basic physics is shown in Fig. 4.2. As a function of the number of molecules 
n in an aggregated cluster, the local free energy /(n) of such a cluster first decreases 
slowly and then faster before increasing again beyond a characteristic size n* . (We 
define /(l) = as the baseline for local packing energy; by 'local free energy' we 
meen a free energy that excludes the translational entropy of the micelle.) The initial 
downward curvature is because although two molecules can lower the energy of their 
hydrophobic tails somewhat by coming together, it is much more efficient to have a 
quorum of molecules so that the tails are entirely separated from water by a layer of 
heads. The upward curvature at n » n* is essentially because for larger n there either 
has to be a hole in the centre of the micelle or it has to become aspherical. If the 
spherical packing is prefered, f(n) is as sketched, but the effect on [i s is no different 
if cylindrical micelles form instead. 

We account for micellization by treating the system as an ideal solution of aggre- 
gated objects (n-mers). According to the Gibbs-Boltzmann distribution the concen- 
tration of each sized aggregate then obeys 

c(n) = v? 1 exp[-/?(/(n) - n/j, s )} (4.9) 

As shown in Fig. 4.2, for strongly negative chemical potential (which holds at low 
enough c s ) surfactant exists primarily as isolated molecules. As n s is raised, c(l) 
increases monotonically. However, a point is eventually reached at which f(n*) — 
n*/j, s = fi s so that c(n*) = c(l). By this point, since n* 3> 1, almost all surfactant has 
formed micelles, and the total concentration c s is about n* times larger than c(l). This 
situation is initially puzzling but has a simple explanation. For [i s <C ^* (where [i* s is 
defined roughly by the condition that the line n* s n is tangent to the f(n) curve) one 
has mostly unaggregated monomers with some exponentially rare micellar aggregates. 
As [i s approaches /i* (so that fj, s = ^* — 6, say), the monomer concentration c(l) 
becomes essentially stuck at 

c(l) = v T x exp[^*] = c* (4.10) 

while any excess surfactant molecules beyond this ceiling form micelles (of size n ~ n* ) 
instead. Because of the form of (4.9), the overall surfactant concentration can be raised 



Micellization 21 




Fig. 4.2 The physics of micellization. (a) Upper curve: the local free energy f(n) as a function 
of aggregation number n has negative curvature at small n and positive at larger n, reflecting 
the existence of a preferred micellar size. The chemical potential /x s is negative but rising 
with c: the three straight lines are examples of the function fj, s n. The concentration c(n) of 
n-mers is the exponential of the vertical separation between /(n) and such a a straight line. 
As Us increases, c(n*) is initially negligible but as Li a — » fx* it abruptly overtakes the monomer 
concentration c(l). (b) The resulting dependence on c s of both c(l) and c mice ;; es ~ c s — c(l). 
(c) The resulting dependence of li s on c s , showing near-saturation at c > c*. 

by a large factor, of order n* , for a very modest (order fc^T) increase in chemical 
potential. Thus, for practical purposes, it is a good approximation to say that /i s 
hits a firm ceiling just below /j,*, beyond which it effectively ceases to depend on c s . 
Likewise the monomer concentration c(l) saturates at c* , which is commonly known 
as the critical micelle concentration or CMC. These behaviours are sketched in Fig. 4. 2. 
In practice there is a rounded corner on jj, s and c(l) rather than a singularity, and 
it is important to note that micellization represents merely a sharp crossover not a 
true phase transition. However, the curvature of these features diverges as n* — > oo. 
In that limit the mathematics of micellization coincides precisely with that of Bosc 
condensation for ideal quantum fluids, as addressed in standard statistical physics 
textbooks [10]. 
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4.2.1 Evolution of the interfacial tension 

According to the above arguments we now have two general classes of behaviour de- 
pending on whether the critical micelle concentration c* lies below or above c. Recall 
that this is the concentration of an ideal surfactant solution at which 7 would effec- 
tively vanish. In case 1, c* < c, so that 7(c s ) follows the ideal curve (Fig. 4.1) so long 
as c s < c* but then abruptly stops decreasing as micellization intervenes: the chemical 
potential fj, s then saturates and no further decrease in tension is possible. In case 2, 
c* > c so that 7 hits zero before micelles are formed. At this point, if water and oil 
are both present in bulk quantities, the system can minimize its free energy by creat- 
ing a macroscopic amount of interface on which the surfactant can reside in comfort. 
When this happens fj, s again saturates: adding further surfactant simply creates more 
surface at fixed fi s . (Hence 7 can never actually become negative as the ideal solu- 
tion calculation might suggest.) By this reasoning, micelles never form under 'case 2' 
conditions. 

The above represents an oversimplified, but useful, picture of the effects of sur- 
factant on interfacial tension. In practice, contributions from curvature energy and 
entropy (considered below) may mean that interface forms spontaneously while 7 re- 
mains slightly positive. Indeed, we defined 7 as the tension of a flat interface; but 
the proliferation of interface occurs when the free energy cost vanishes for creating it 
not in that state, but the state of mimimum free energy (which might be crumpled 
or curved). Nonetheless, it is broadly correct to distinguish between case 1, where 7 
is reduced considerably but remains or order its original value, and case 2 where 7 
becomes effectively, if not actually, zero. (It is probably not correct to assume that 
micelles never form in case 2, however.) For a far fuller molecular discussion of how 
surfactants modify interfacial tension, see [26]. 

4.3 Finite tension: metastable emulsions 

Because c* is generally small (typically 10~ 2 Mol 1 _1 , often far less), case 1 is generally 
the more common: the effect of surfactant is to reduce interfacial tension to half or a 
third of its previous value. Because 7 remains finite, the global minimum of free energy 
is always that of coexisting A-rich and B-rich phases, separated by a flat interface. The 
area of this interface is set by the container geometry (and, in practice, gravity) and 
scales as V 2 ^ 3 where V is the sample volume. Interfacial physics thus contributes 
negligibly to equilibrium states in the thermodynamic limit. Nonetheless, one can 
create metastable emulsions, for instance by stirring. These are generally droplets (of 
A in B, say) but by drainage under gravity, for instance in a centrifuge, much of the 
continuous B phase can often be expelled to create a so-called biliquid foam [3]. 

As the name suggests, biliquid foams are very similar in structure to foams made of 
air bubbles in a surfactant solution (soap froths). They comprise polyhedral droplets 
of A (say) separated by thin films of B, and the stability of the foam depends on a 
barrier to coalescence of A across these thin films. The role of surfactant is more to do 
with raising this barrier (via coulombic, entropic or steric forces in some combination) 
than with reducing 7. In many cases, biliquid foams can persist for hours or days, and 
sometimes longer. To achieve this one must suppress not only the rupture of thin films 
but also the Ostwald process which, despite the more complicated geometry, still drives 
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diffusion of A from small (few-sided) to large (many-sided) polyhedral droplets [27; 
28] . Inclusion of a trapped species helps, as described previously, but this must now be 
extremely insoluble in B so as to have negligible diffusion even across the thin B films 
present in the foam structure. So long as they remain metastable against rupture and 
coarsening, biliquid foams, like soap froths, are solid materials (generally amorphous, 
though ordered examples can be made). As such they have an elastic modulus, and 
also a yield stress, both of which scale as G <~ j/R with R the mean droplet size. 
This is an interesting example of a solid behaviour emerging solely from the spatial 
organization of locally fluid components - for even the surfactant on the interface is 
(normally) a fluid film. 

4.4 Effectively zero tension: stable microemulsions 

We return now to to case 2, where a sufficient level of added surfactant can reduce 
7 to negligible levels for c s > c. This can lead to thermodynamically stable emul- 
sions, generally called "microemulsions" . As previously described, once this happens, 
enough A-B interface is created to accommodate all surplus surfactant, of which the 
concentration is c s — c. Since c lies below the CMC, which is generally itself small, the 
interesting range is usually c s > c so that one can treat effectively all the surfactant 
as interfacial. The interfacial area S of the fluid film then obeys 

^ = (c s -S)E~c s E=^E (4.11) 

Here E is a preferred area per surfactant molecule; <fr s is the volume fraction of sur- 
factant and v s its molecular volume. Note that principle, the area per molecule could 
deviate from its preferred value - an effect important when studying, for instance, 
Langmuir-Blodgctt films of insoluble amphiphiles such as lipids. However, the soluble 
surfactants normally used for emulsification have binding energies to the interface that 
are only modest (say 5 — 12/csT), so they can adsorb and desorb from the interface on 
a short timescale. These processes can rapidly restore the preferred area per molecule. 

Treating E as constant, the specific interfacial area S/V is fixed directly by <p s 
via (4.11). We next ask, what is the configuration of the interface? This is set by 
a competition between entropy, which prefers small, wiggly structures, and bending 
energy, which prefers extended, smooth ones. The bending energy can be treated 
by a leading order harmonic expansion about a state of preferred curvature that is 
set by the molecular geometry of the surfactant layer (and is tunable by varying that 
geometry, or by mixing surfactants with different curvature preferences). By a theorem 
of differential geometry, at each point on the A-B interface one can uniquely define two 
principal radii of curvature, R\,R2- These are the radii of two circular arcs lying in 
perpendicular planes to one another, each of whose centres lie on the surface normal, 
such that each arc kisses the surface without crossing it; see Fig. 4. 3. (Alternatively, 
Ci : 2 = 1/Ri,2 are the eigenvalues of the surface curvature tensor; see F. David's 
chapter in [5].) The radii are signed quantities and we shall take them positive for 
curvature towards A. For a spherical droplet of A of radius R, we have R\ — Ri = R 
whereas for a cylinder of radius R, we have R\ = R and R2 — 00. A saddle shape has 
i?iand R 2 of opposite signs. 
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Fig. 4.3 The grey lines are intended to represent a curved surface such as a piece of bicycle 
tyre. At the point whose surface normal is indicated by n, the construction of the principle 
radii of curvature Ri and R2 is shown. 

The harmonic bending energy then reads 




(4.12) 



There arc 3 material parameters, K, K and Rq. Note that a general expansion to 
second order in two curvatures would have five parameters {a,\C\ + a 2 Cf + a 3 CiC 2 + 
(X4C1 + 0502) but for a fluid film a\ = 0,2 and (Z4 = as by rotational symmetry; the 
remaining expression can be reorganized to give (4.12) [5]. The chosen parameters K 
and K have dimensions of energy while i?o is a length defining the preferred radius of 
mean curvature. Theories exist to relate these quantities to the molecular geometry of 
surfactants [4], but we do not pursue these here. 

We now note that the interface between A and B cannot end in mid air: any edges 
must lie at the boundary of the container. (With periodic boundary conditions, no 
edges are possible.) The interface S can have disconnected parts (droplets) but must 
be orientable so that A is enclosed by it and B excluded. Moreover, the volume Vi n 
enclosed by the interface obeys 

^=d> A + | = $ (4.13) 

Here we have partitioned the surfactant equally between A and B to allow us to define 
the volume Vj n as enclosed by a mathematical surface of no thickness; the phase volume 
of Vin is then $ (with V ou t = 1 — $), and a completely symmetric state has $ = 1/2. 

To examine the statistics of the interface and determine its free energy, we now 
ought to compute F = —k B T\nZ where the partition function 
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Z = / cxp[-f3F bend }VS 



(4.14) 



is found by integrating the Boltzmann weight computed using (4.12) over all surfaces S 
that enclose volume $>V. This has been an active area of statistical mechanics for over 
30 years, and we have time here only for a brief flavour of the topic. Approaches range 
from crude estimates (exemplified below) to what is essentially string field theory. 
Although our problem is not quantum mechanical, it shares with string theory the 
thorny issue of how exactly to count the continuum of distinct configurations accessible 
to a two dimensional fluid manifold embedded in a higher dimensional space (here) or 
space-time (strings) (see F. David in [5]). 

4.5 Some useful concepts relating to bending energy 

Here we summarize three useful concepts, using which large parts of the statistical 
physics problem referred to above can be qualitatively understood. 

4.5.1 Gauss-Bonnet theorem 

This theorem states that 



here N c is the number of components of our surface (where a component is a dis- 
connected piece such as a droplet) and Nh is the number of handles. A handle is a 
doughnut-like connection between one part of the surface and another. Thus for a 
sphere N c = 1 and N h = whereas for a torus, N c = 1 and Nh = 1. Accordingly the 
bending energy term governed by K in (4.12) vanishes for a torus but not a sphere. 
Importantly, the result is topologically invariant so that any deformation of a torus 
still has zero for this quantity and any closed droplet that is deformable continuously 
into a sphere has the same value of 4irK . Thus the bending constant K contributes a 
term to F bend that does not care about the local deformations of the surface, only its 
topology. 

In surfactants that stabilize emulsification, K is generally negative. To understand 
this, one needs to be aware of the existence of periodic surfaces of constant mean 
curvature. These comprise a periodic surface element (Fig. 4. 4) which connects with 
identical copies of itself in neighbouring unit cells to create a structure with only one 
global component, but several handles per unit cell. Such surfaces can be made entirely 
of saddles having the required mean curvature i?o (at least when R is large so the 
preferred curvature is weak or negligible) . Accordingly the K term in the bending free 
energy (4.12) vanishes everywhere. Because the K term vanishes, all depends on K. 
If this is positive, handles are favoured, and this means that the periodic structure 
would like as small a unit cell as possible. (Its shrinkage is ultimately controlled by 
anharmonic terms in the free energy that were not included in (4.12).) The result is a 
short length-scale triply periodic liquid crystal, which is quite an interesting structure 
in itself [29] but not an emulsion as such. Thus we assume K < in what follows. 
We shall also assume 2K + K > so that the bending energy of a sphere is positive 
(otherwise one expects instead a proliferation of tiny spheres). 




(4.15) 
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Fig. 4.4 A sphere, a torus, and the unit cell of a periodic surface of constant (approximately 
zero) mean curvature. The hole through the torus is a handle. The grey discs on the periodic 
surface are cuts across it at the junction points between unit cells. Gluing a pair of these discs 
together at the faces of the unit cell creates one handle. Thus the final periodic structure has 
three handles per unit cell, but only one global component (since the entire surface becomes 
a single connected object). 

4.5.2 Persistence length 

We now set i?o = oo so a flat interface is preferred. The bending energy can then be 
evaluated for small fluctuations in shape described by a height field h(x, y) above a 
flat reference plane. One finds [30] 

F bend = J (^(V 2 h) 2 ^dxdy~*J2q 4 \h q \ 2 (4.16) 

where in the first expression V 2 is defined with respect to the x and y coordinates and 
in the second we have taken a fourier transform of the height field. 
Equipartition of energy then demands 

<N 2 > cc ^ (4.17) 

which holds for q < ir/£ with I some cutoff length comparable to the thickness of the 
surfactant film. From this it is a simple exercise to show that 

<|V%)-VM0)| 2 )cxM!ln(£) (4 . 18) 

Thus the orientation of the surface deviates from its initial value with a logarithmic 
dependence on separation; when this deviation is large, the expansion underlying (4.16) 
breaks down. An equivalent statement is to introduce a persistence length [31] 



£,k — £ exp 



4ttK 
akgT 



(4.19) 



where a is sometimes referred to as a 'geometrical' constant. The interpretation is that 
a fluid interface governed by (4.12) completely loses its sense of orientation at length 
scales beyond £k ■ 

Because of the exponential factor, £k is exquisitely dependent on K (and also on a, 
which turns out to obey a — 3 as discussed below). Thus for I = 1 nm, £ K ~ 1/xm when 
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K = 1.65/csT. For K/ksT = 3, we already have > 300^m, and is irrelevantly 
large, for our purposes, once KjksT is much larger than this. (Indeed for K/ksT = 6 
we have £k = 9cm.) 

4.5.3 Renormalization of the bending constant 

In many cases we do not want to describe the A-B interface in complete detail, but 
instead need a coarse grained description of its properties on some geometric length 
scale £ set by, for instance, the size of droplets in an emulsified state. Under coarse 
graining we replace an entropically wiggly interface, with structure at short scales 
characterized by the bending spectrum (4.17), by a smooth one on the scale £ for 
which such structural detail is still present but removed from the description. To allow 
for that, we need to accept that the effective bending constant for the smoothed-out 
interface is softer than the microscopic one. Indeed the result of a careful and quite 
technical field-theory calculation (see F. David in [5]) is that, to leading order, 

tf e/ ,=*-^]n(f) (4.20) 

This is quite plausible in that it suggests there is effectively no resistance to macro- 
scopic curvature at length scales beyond £k as defined in (4.19). Indeed consistency 
with this calculation requires the choice a = 3 as previously mentioned. However, the 
numerical prefactor in (4.20), which indeed fixes this value of a, is far from obvious 
and was the subject of heated debate for a number of years. 



4.6 Some consequences of bending energy physics 

In addressing this large topic, we will have time for only a few limiting cases. 
4.6.1 Emulsiflcation failure 

We first consider K > 0, and I -C Rq <C £k- This holds for interfacial films with finite 
preferred radius of curvature that are relatively stiff on that length scale. Entropy is 
negligible: we need only minimize Fbend at fixed SV . The computation of Fb en d for 
spheres, cylinders and lamellae is trivial and for simplicity we limit attention only to 
these geometries. We have 



Fbend = 4ir[2K 
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(4.21) 



for a sphere of radius R; 
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(4.22) 



for a cylinder of radius R and length L; and Fbend = 2KAjR\ for a flat sheets of 
area A. A schematic phase diagram found by comparing these forms [4] is shown in 
Fig. 4. 5. This shows, as one might expect, a preference for lamellae when R is large 
and spheres when it is small. More interestingly, it is a simple exercise to show (by 
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equating the enclosed volume to V<fr and the surface area to S) that the droplet size 
R for (monodisperse) spheres obeys 




(4.23) 



where I — v s /T, is the thickness of the surfactant film. (Here v s is an effective molecular 
volume for surfactant.) If R s < Ro, then, modulo slight shifts relating to the value of 
K, the droplet phase is stable. However, if surfactant is removed or the internal phase 
volume fraction $ is increased to the point where R s found above exceeds Ro, there is 
no advantage to the system in paying the additional bending cost of having droplets 
larger than their preferred curvature radius. Instead, the droplets remain of size R , 
so that the value of <E> in the droplet phase is less than its globally imposed value. This 
means that an excess phase comprising bulk A-rich fluid is expelled from the system. 
The system has created the perfect interface (in area and curvature) for the amount 
of surfactant present, but this fails to enclose all of the dispersed phase - a situation 
known as emulsification failure [4; 32]. 




-K/ (XK ^ 

Fig. 4.5 Schematic phase diagram (after [4; 32]) showing the phase of lowest free energy as a 
function of enclosed phase volume the elasticity parameters Ro, K, K, and the film thick- 
ness t. This is a simplified calculation which does not account for entropy nor the existence 
of miscibility gaps between the states shown. Nor does it include anharmonic corrections to 
the bending energy. 



4.6.2 Bicontinuous microemulsions 

The second case we address here is where Rq for simplicity we treat i?o as 

infinite so the preferred interfacial state is symmetric between A and B (sometimes 
called a 'balanced' system). We also assume £ -C £k -C 100/im, so that entropy matters 
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(including the renormalization of K) and can compete with bending energy on roughly 
equal terms. 

Assuming <£> of order 0.5 (roughly symmetric amounts of A and B) we can introduce 
a structural length scale £ which is then set by </>„. Specifically for a lamellar phase 
one has a layer spacing £ between adjacent surfactant films set by 

S - ^ (4.24) 

(This is subject to logarithmic corrections from fluctuations which need not concern 
us here.) When £ ~ I the system has no option but to fill space with flat parallel 
layers. As cf) s is reduced (£ raised) the layer spacing £ becomes comparable to £k- 
For £/£,k < 1/3 (or so), the lamellar phase fluctuates but remains stable. However 
if 4> s is then decreased further so that £ ~ these layers melt into an isotropic 
phase comprising (for <j> ~ 0.5) bicontinuous domains of A and B fluids separated by a 
fluctuating surfactant film (Fig. 4. 6) [31]. This is the bicontinuous microemulsion and 
represents a thermodynamic route to prevent coarsening of the transient bicontinuous 
emulsion structures encountered in Section 3. If $ now deviates strongly from 0.5, then 
(just as found there) the structure depercolates, forming a droplet phase; the details 
of where and how this happens involve K which, we recall, is sensitive to topological 
changes. The droplet phase so formed is somewhat different from the one discussed 
above at Rq -C £,k, since this one is stabilized by entropy and fluctuations, not by a 
preferred curvature of the droplets. 

Several theories of the bicontinous microemulsion were developed in the 1980s. 
Some of the most successful [33; 34] used coarse grained lattice models in which fluid 
domains are placed at random on a lattice of some scale £ ; the bending energy and area 
of the resulting interface can be estimated and used to calculate a phase diagram. One 
specific feature is the appearance of three-phase coexistence in which a 'middle phase' 
microemulsion coexists with excess phases of both oil and water. This is somewhat 
analogous to a double-sided emulsification failure in which, for a fixed quantity of 
surfactant, the system creates a happy interfacial structure on a certain length scale 
and then rejects any excess oil and/or water. This length scale is now set by rather 
than Rq. Although the models do predict this, it remains somewhat unclear even now 
why this structure does not want to fragment further. Such fragmentation would lead 
to two coexisting phases of dilute A droplets in B and vice versa. One possibility is 
that K e ff(£), whose renormalization properties we have not discussed (see F. David 
in [5]), becomes positive at £ > £k- This would cause condensation into a handle-rich 
state, similar to the periodic ones previously stated to arise for K > 0, but stabilized 
now by an entropy-driven tendency to form handles that only operates at large enough 
length scales. 

4.6.3 The sponge phase 

One final twist on this story concerns systems in which there is a huge phase volume 
aysmmetry between A-rich and B-rich fluids, but where the surfactant has a very 
strong molecular preference to form a flat film rather than highly curved structures 
such as micelles. The interfacial structure that forms spontaneously at 7 ~ is, in 
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Fig. 4.6 Schematic evolution from lamellae to bicontinuous microemulsion as the concen- 
tration of surfactant (and hence £,k /£) is decreased in a system with * ~ 0.5 (three central 
frames). The breakdown into droplet phases occurs on varying * away from 0.5, with details 
dependent on K. 

the almost complete absence of B, necessarily now a bilayer with A (usually water) on 
both sides and a thin B layer in the middle. It is easy to imagine the lamellar state, 
and for small volume fractions of bilayer this suffers the same instability towards 
melting as described above for the microemulsion. The result is a bilayer film that 
now separates two randomly interpenetrating domains containing the same solvent A: 
this is usually called the 'sponge phase'. Perhaps surprisingly, the quantity $ which 
discriminates between Vi n and V ou t remains meaningful. However it is now fixed not 
by the global phase volumes of the two solvent domains (since these contain the same 
solvent and are interchangeable) but instead takes whatever value minimizes the free 
energy. Moreover, defining = 2$ — 1, the microscopic free energy of this system is 
invariant under ^ o — * (an operation which exchanges the identities of the 'inside' 
and 'outside' fluid domains). This symmetry can remain intact in the physical state of 
the system, or it can break spontaneously, for instance by the formation of a phase of 
closed spherical vesicles [35; 36]. The symmetry breaking can be continuous, and then 
shows critical behaviour of Ising-like character and a unique scattering signature [37] . 

This scenario is theoretically quite rich and elegant [38], but the sponge phase 
contains at most a few percent of its second solvent (B, as described above) sequestered 
within the bilayers. It is therefore not particularly useful for cmulsification and has 
found relatively few applications in technology so far. This contrasts somewhat with 
our next topic. 
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A typical soluble surfactant has a hydrocarbon tail containing 10-20 repeat units 
and a head group with either a dissociable salt (ionic or anionic surfactants) or a 
polar polymer such as a string of ethoxy groups (typically 4-8 of these in a nonionic 
surfactant) . Such surfactant molecules generally have energies of attachment to the oil- 
water interface of between 5ksT and 20/cb? 1 , high enough to alter interfacial properties 
but low enough to allow rapid adsorption and desorption to maintain local equilibrium. 
However, it is easy to create amphiphilic species that are much larger [40; 41]; these 
range from lipids, via block copolymers (tens to tens of thousands of repeat units) and 
globular proteins, to so-called 'janus beads'. The latter are colloidal spheres, up to a 
micron in size, with surface chemistry that favours water on one hemisphere and oil on 
the other. For typical solid-fluid interfacial tensions (~ 0.01 Nm~ 2 ) janus beads have 
attachment energies of order 10 7 fcsT or larger. Such species are adsorbed irreversibly in 
the sense that Brownian motion will never lead to detachment. All of these amphiphilic 
particles can, in principle, lead to a vanishing of the thermodynamic interfacial tension 
if micellization and similar processes are avoided. However, thermodynamic concepts 
such as interfacial tension itself turn out to be of very limited value in understanding 
irreversibly adsorbed layers. 

5.1 Adsorption of non-amphiphilic colloids 

In fact, although janus particles can be made and arc often studied [39], they are 
rarely used in practical applications to stabilize emulsions. This is because a simpler 
route exists, in which these purpose-built amphiphilic colloids arc replaced by those 
of homogeneous surface chemistry: that is, ordinary colloidal spheres [6; 42; 43]. 

Such particles can also have attachment energies to the interface that are vastly 
larger than ksT. This is most easily seen when the surface chemistry is chosen with 
equal affinity towards water or oil so that the two solid-fluid interfacial tensions, 750 
and "fsw are the same. The energy of such a particle (of radius a) is independent of 
where it resides in an oil-water system, but the energy of the oil water interface is 
reduced by ira 2 -f if the particle is placed there. This is simply the surface energy of the 
interfacial patch now covered by the colloid, and 7 is the usual fluid-fluid interfacial 
tension. Setting 7 = 0.01Nm~ 2 gives AF attach ~ 10 7 fcsT for a = 1/xm and AF attach ~ 
lOksT for a = lnm. Thermal detachment thus remains negligible for a > 5nm. 

More generally, the two solid-fluid tensions are different. However, the lowest free 
energy state has the colloid on the interface so long as the contact angle 6, which obeys 



7 cos6> = 7s - -ysw 



(5.1) 
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obeys < 9 < ir. (This is the definition of partial, as opposed to complete, wetting.) 
In the absence of a body force such as gravity acting on the particle, the interface 
remains perfectly flat; the particle is displaced so that it intersects the interface at 
angle (Fig. 5.1). Accordingly the area of the covered disc is reduced; the attachment 
energy is also, but remains very large compared to fceT, at least for a > lOnm, unless 
cos 9 is very close to ±1. Indeed, a calculation of the force Jd required to detach the 
particle gives [6] 

f D =7ra 7 (l±|cos0|) (5.2) 

with the +(— ) sign applies when the particle is being pulled out of (into) its preferred 
solvent. Modest deviations from the neutral wetting condition (co&9 = 0) thus have 
little effect on the detachment barrier. 




Fig. 5.1 Geometry of a partially wet colloidal particle at the fluid fluid interface; 9 is the 
contact angle. This is usually measured through the polar phase, so the upper fluid is water 
and the lower oil as drawn here. 

Turning now to the case of a layer of adsorbed particles, it is often said that devia- 
tions from neutral wetting creates spontaneous curvature, as was described previously 
for surfactants. (In principle this effect looks very strong, with preferred curvature 
radii of order a/cos9.) However, this is deceptive, since the particles are perfectly 
happy on a flat interface unless their density is high enough to jam them into contact. 
Beyond this point, it is true that the interfacial area can be reduced further at fixed 
particle coverage by introducing curvature, but without any constraint on the overall 
interfacial geometry it is even better to expel the particles and have no interface at 
all. Usually the overall geometry is constrained purely by the necessity of having an 
interface between bulk coexisting phases; in that case, regardless of 9, the lowest free 
energy state is found first by creating a flat interface, and then by covering as much 
of it as possible with particles. At no point is a curved interface actually preferred. 

5.2 Particles on curved interfaces 

Consider now placing a number of colloidal particles on a spherical droplet of fluid A 
in fluid B, or vice versa. It is crucial to note that each colloid can be accommodated 
with the required contact angle 9 simply by cutting a small spherical cap out of the 
interface and slotting the particle into place there. The contact line is a perfect circle, 
as required for tangency at fixed angle to a sphere. To conserve the enclosed volume 
the droplet radius may change slightly, but it remains perfectly spherical. There is no 
change in Laplace pressure, and the energy is completely independent of where the 
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colloids are placed. Accordingly there is no (tension-mediated) capillary force between 
the particles. These statements can only change if particles become jammed so they 
interact directly via particle-particle forces. 

Contrast this with the case of spherical colloids on a hypothetical cylinder of fluid. 
It is not possible now to insert a spherical particle into the surface of this cylinder at 
fixed contact angle 6 unless the fluid interface becomes deformed. This deformation 
costs extra interfacial energy, and can be minimized by placing two particles close 
together rather than far apart. Accordingly, there is a capillary attraction between the 
colloids, which will have a strong tendency to aggregate. Similar arguments apply to 
nonspherical particles such as ellipsoids, even on flat surfaces [44; 46]. More generally, 
this field offers various formal mathematical problems involving area minimization 
under the constraints of enclosed volume and adsorbed particle number, many of 
which seemingly remain unexplored. 

5.3 Particle-stabilized emulsions 

The use of partially wetting colloidal particles to stabilize emulsions dates back at 
least one century [42]. Unless janus particles are used [45], the resulting 'Pickering 
emulsions' are always metastable: as stated already, the minimum free energy state 
comprises bulk A and B phases, separated by a flat interface, with as much of this 
interface as possible covered by particles. (This requires an ordered hexagonal packing 
in principle.) The remaining particles are then distributed randomly in their preferred 
solvent (if 9 7^ 7r/2) or in both solvents if there is no preference (6 = it/2). 

However if the system is stirred or otherwise agitated, it is quite easy to create 
robust droplets, whose size is fixed by the production process as well as the fluid 
phase volumes and the particle density. Closely related droplet phases can also be 
made in which the dispersed phase comprises air or vapour rather than liquid [47] . 

5.3.1 Resistance to coalescence 

One specific route to Pickering emulsions [6; 48] is to create a vast number of very 
small droplets by applying extreme flow conditions to temporarily mix fluids A and 
B. The flow sweeps particles onto the interface whose initial surface area is, however, 
much larger than they can cover. Coalescence initially proceeds as normal (assisted by 
maintaining a lower but finite flow rate) . During a coalescence between two droplets of 
size R, the total volume they enclose is conserved, so the radius changes as R — > 2 1 / 3 i? 
and the surface area as S — > 2~ 1 / 3 <S. Thus the area fraction of particles increases 
by 2 1 / 3 each coalescence until the surface density is high enough to prevent further 
coalescence. 

Because it requires significant local rearrangement of the particles in relation to 
the interface, the resulting energy barrier scales as na 2 -y ^> k B T. (Indeed, Pickering 
formulations often offer exceptional stability against coalescence compared to conven- 
tional surfactants.) This requires a coverage of comparable to, but significantly below, 
that of a densely packed 2D amorphous film. 

In most cases the protected droplets therefore have a fluid layer of particles at 
their surface not a jammed one, so that they relax to a spherical shape. In some cases, 
though, two droplets coalesce that are already nearly protected. The reduced surface 
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area of the combined droplet, if spherical, is then not enough to accommodate all the 
particles. Coalescence proceeds, but then is arrested in a jammed state comprising a 
non-spherical droplet, whose surface particles are clamped in position by interfacial 
tension [49]. (Related effects can be seen at lower particle coverage if there are strong 
attractive forces between colloids creating a bonded rather than a jammed interfacial 
layer.) Similar jammed structures can again been seen in armoured air bubbles [47]. 

5.3.2 Resistance to Ostwald ripening 

Recall that a close packed monolayer of particles (whether ordered or amorphous) can 
be placed on the surface of a spherical fluid droplet of radius R (say) without altering 
its interfacial geometry. Imagine such a droplet with the particles just in contact with 
one another. There is still a fluid-fluid interface at the interstices between particles, 
and this has the same curvature, and hence Laplace pressure j/R, as the original drop. 

Suppose now that this droplet is in diffusive equilibrium with one or more larger 
ones. According to the Ostwald mechanism, it will start to shrink. However, if the 
particles are already in contact they cannot follow the droplet surface inwards as 
this happens. Moreover, each particle demands an unchanged contact angle with the 
interface, which is effectively now pinned to the particle layer. It is easy to see (Fig. 
5.2) that even a small loss of volume of the droplet under these conditions will cause an 
elimination of the Laplace pressure followed by its sign reversal. Thus a state of finite 
droplet radius, but zero mean curvature of the interface so that AP = 0, presumably 
exists. (Proof of this is one of the seemingly unsolved problems referred to in Section 
5.2 above.) The resulting droplet is fully resistant to Ostwald ripening, and therefore 
a suspension of Pickering droplets is also stable against it, so long all the interfacial 
particle layers are jammed. The mechanism is similar to that described in Section 3.6 
using a trapped species, with the interfacial particles themselves playing that role. 

If the jamming is not complete initially, as would typically be true for emulsions 
made by the arrested coalescence method outlined above, a minority of droplets may 
grow, ultimately forming a coexisting bulk phase of the dispersed fluid. However, if 
this is now removed, the emulsion component is highly stable against both coalescence 
and Ostwald ripening. There arc subtleties with this simple picture, though, if the 
dispersed phase contains more than one soluble component [50]. 

5.4 Particle detachment 

The barrier to individual particle detachment ceases to be very large compared to fc^T 
for nanometre scale particles. (With such particles, the interfacial energies we have 
described may need supplementing with an energy term proportional to the length of 
the contact line [43].) The same can apply for larger colloids if the interfacial tension 
7 is relatively small, or if surfactants are added that change the particles' wettability. 
Low tension often arises for polymer polymer mixtures as opposed to oil and water; 
and blending of polymers into emulsions is an important area of technology. Moreover, 
once particles are jammed by interfacial tension into a dense layer, the detachment 
is a collective process and it is not clear whether some pathways then exist that 
allow particles to escape without crossing a high barrier [51]. Addressing the physics 
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Fig. 5.2 In light grey is the initial locus of a fluid interfaces into which are inserted the 
particles shown in darker grey. These are jammed in a 2D layer but have interstitial fluid 
regions as shown between them. If the volume of fluid in the droplet is now reduced, to 
maintain a fixed contact angle with particles that cannot move, the curvature of the interface 
is reversed to give the final fluid locus (black) . This creates a negative Laplace pressure which 
switches off the Ostwald process. 

of Pickering emulsion with non-negligible detachment processes is a complicated but 
interesting open area. 

5.5 Some interesting particle-stabilized structures 

More generally, there are large gaps in our theoretical understanding of particle- 
stabilized emulsion formulations, in part due to the fact that non-desorption causes 
almost all the interesting structures to be far from equilibrium. In what follows we 
describe a selection of these structures, as depicted in Fig. 5. 3. In at least one of these 
cases, simulations have been a guide to the creation of new states [52]. 

5.5.1 Pickering foams 

Drainage or ccntrifugation of a Pickering emulsion can lead to a compressed foam 
structure that is like the biliquid foams discussed previously [53]. These can be quite 
stable thanks to the combined resistance to coalescence and Ostwald ripening pro- 
vided by the particles. However, removal of the continuous phase causes droplets to 
deform into polyhedral shapes which increases their area and can bring the particle 
coverage below the threshold for protection. Interestingly, for cos 9 ^ it is possible 
for a particle to bridge across a thin film of its preferred solvent while maintaining 
the equilibrium contact angle at both interfaces. Though not as effective as two mono- 
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Fig. 5.3 Various particle-stabilized interfacial structures between fluids A and B. Notation: 
the continuous squiggly line stands for a monolayer of colloidal particles, (a): A biliquid foam 
with thin B films separating polyhedral A droplets. Inset: the possibility of 'zipping', (b) A 
multiple emulsion, (c) A bijel; compare the bicontinuous microemulsion in Fig. 4.6. 

layers, the resulting 'zipped interfaces' still present a barrier to film rupture [54]. If 
present, these might also resist re-expansion of a foam after centrifugation has ended. 

5.5.2 Multiple emulsions 

Simple manual agitation of binary immiscible solvents containing partially wettablc 
particles often results in droplet-within-droplet structures known as multiple emul- 
sions. Such structures require stability against both coalescence and ripening for A 
droplets in B and B droplets in A simultaneously. This is relatively difficult for sur- 
factant formulations but seemingly quite easy with particle-stabilized ones. 

5.5.3 Bijels 

Bijels (bicontinuous interfacially jammed emulsion gels) are metastable analogues of 
the bicontinuous microemulsion: a particle layer resides at the interface between inter- 
penetrating domains of A and B. This structure was predicted first computationally 
in 2005 [52] and confirmed in the laboratory in 2007 [55]. It has no known counterpart 
among metastable surfactant-stabilized formulations, presumably because surfactant 
desorbs too easily ever to prevent the coarsening of an interface with finite 7. In a bijel, 
the interfacial film of non-detachable particles is clamped by tension into a 2D jammed 
layer, which imparts solidity to the whole 3D structure [56; 57]. This robustness can 
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be improved further by having an interaction potential between particles with a steep 
barrier and then an attractive minimum at short distances. The interfacial tension 
pushes particles over the barrier creating a permanent interfacial film [58]. 

One recipe for making a bijel is to choose a fluid pair A+B that are misciblc 
at high (or perhaps low) temperature [59]. The colloidal particles are then dispersed 
within the single-phase AB mixture. On quenching the temperature, the fluids separate 
and particles are swept onto the interface. The coarsening process arrests when a 
jammed monolayer is formed, creating the bijel. The final structural domain size obeys 
£ ~ a/<f> p with a,(j) p the particle size and volume fraction, and the elastic modulus 
of the solid bijel scales as G ~ 7/£. Bijels are currently being explored for various 
applications including materials templating [60], where their A-B bicontinuity, and in 
some cases tricontinuity (allowing also for the percolating particle layer) could prove 
advantageous. 

5.6 Stability against gravity 

I choose as a final topic something representative of very recent research done in my 
own group at Edinburgh (with experiments led by Paul Clegg). Consider a droplet of 
oil in water (which generally will float) stabilized by particles of higher mass density 
than either fluid. Each particle feels the force of gravity: for small particles this is 
very weak (though it can be increased by centrifugation) . As we will see, the force 
can accumulate across particles to give much bigger effects than expected for a single 
particle. The same applies when particles feel body forces that are not gravitational 
in origin - for example if emulsion stability is provided by magnetic colloids which are 
then subjected to a field gradient [61]. 

5.6.1 Critical Bond number 

As shown previously in (5.2), a force fo is required to remove one particle of radius a 
from a flat interface. The same force is needed, to a good approximation, for detach- 
ment from a curved droplet of radius R, so long as R >• a. Equating fn to the force of 
gravity f g gives a critical value of Bo*(l) = 1 for the Bond number which (absorbing 
a factor 4/3 often found in the literature) we define as 

4a?gAp 

B ° = 37(1- cos (?) (5 ' 3) 

Here g is gravity and Ap the density mismatch. 

Our task is to generalize this result to find the critical Bond number Bo*(iV) for 
detachment from a droplet coated with a monolayer of N particles. (This will depend 
on the exact coverage, which we assume to be constant and fairly close to a packed 
monolayer.) Because the gravitational force can be transmitted from one particle to 
those below that support it, it is a reasonable expectation that Bo*(iV) < 1. If so, how 
does it depend on ./V? Two possible mechanisms come to mind, as follows. 

5.6.2 Keystone mode versus tectonic mode 

The first candidate is that a finite fraction of all the weight of the particles bears down 
on a 'keystone' particle at or near the bottom of the droplet. The net gravitational 
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force is N times that on one particle, while fo is effectively unchanged. The keystone 
particle detaches once Nf g > fo] this implies that Bo*(iV) cx N -1 . A stream of 
detachments will continue until the number of remaining particles falls below Jd/ fg- 
This can leave the upper part of the droplet unprotected against coalescence, resulting 
in macroscopic instability of the emulsion. 

The second candidate mechanism is that a finite fraction of the droplet breaks off 
in a collective 'plate tectonic' detachment event. A crude estimate of when this might 
happen is found by balancing the gravitational force on the lower half of the droplet 
Nfg/2, against the interfacial tension across the equator, 2ttjR. Since R ~ aN 1 ' 2 this 
gives (within prefactors that depend on geometric details such as the deformed droplet 
shape) Nf g ~ N 1 / 2 f D , so that Bo*{N) cx AT" 1 / 2 . 

For droplets containing tens to thousands of particles, the competition between 
these two modes seems to be finely balanced. Experiments on particles with signifi- 
cantly non-neutral contact angles 9 suggest the keystone mode [62] whereas simulations 
with 9 = w/2 suggest the tectonic one [63]. So a provisional conclusion is that the re- 
sult depends on contact angle; this is plausible because the one-particle detachment 
energy (at given particle size) depends strongly on contact angle, whereas the force 
balance in the tectonic mode does not. 
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